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The dynamics of the interface between two immiscible fluids in a rotating Hele-Shaw cell are stud- 
ied experimentally, theoretically and by phase-field simulations of the H-S equations using standard 
boundary conditions. As the central, denser fluid is centrifuged, it forms fingering patterns with 
long, thin radial filaments ended by a droplet, alternating with incoming fingers of the outer, less 
dense fluid. Simulations show the length (width) of the filaments to grow (decay) roughly exponen- 
tially, and the incoming finger tips to asymptotically approach a finite radius for n-fold symmetric 
initial conditions; these thus tend to a stationary-shape "star-fish", whose form is calculated. The 
filament width decays with a time constant which depends only on the viscosity contrast, whereas 
its length exhibits a completely universal growth rate; the latter is related to the run away of an 
isolated droplet, for which we give an exact solution. The exponential behavior is clear for high 
viscosity contrasts A, but important deviations are found for low A. Both experiments and simula- 
tions show systematic pinch-off of the droplets at the tips of the filaments for low and not for high 
A. A lubrication approximation is derived and successfully accounts for the filament thinning and 
the differences with A; in particular, it explains why pure exponential thinning is not observed for 
low A, and it could clarify the presence or absence of finite-time pinch-off itself, since the agree- 
ment of experiments and simulations suggests that this phenomenon is contained in the Hele-Shaw 
equations. This agreement includes both high- and low-A morphologies, and the growth rate of 
the filament length for high A; for low A, the experimental time constant appears to be different 
from that predicted by standard Hele-Shaw boundary conditions and observed in simulations. An 
effective slip condition for the Poiseuille flow of inner liquid across the cell gap in the case of two 
liquids gives a possible explanation of this discrepancy. 

PACS numbers: 47.55.Dz, 47.20.Ma, 47.ll.-fj, 68.03.Cd 



I. INTRODUCTION 

Topological singularities such as interface pinch-off in 
fluid flows have been the object of intense study in the 
last decades 0,1111 0,1110,11113111 13 13 111 The 

pinch-off of a droplet from a volume of fluid embedded in 
another one is easy to observe in everyday life, and has 
fascinated both physicists and applied mathematicians. 
Because the interface between the two fluids shrinks to 
arbitrarily small scales at the point of detachment, a con- 
tinuum, coarse-grained description must break down at 
some stage. The failure of the macroscopic description is 
reflected in the spontaneous generation of singular behav- 
ior. In the neighborhood of such singularities, the prob- 
lem might become scale-free, giving rise to self-similar 
scaling behavior, for which some degree of universality 
is naturally expected 0. Interestingly, hydrodynamics 
alone have been found sufficient to account for the oc- 
currence of finite-time pinch-off in several cases. It is 
remarkable that the continuum, hydrodynamic descrip- 
tion may correctly predict the occurrence of finite-time 
singularities out of a smooth initial condition, and prop- 
erly describe them all the way to the very microscopic 
scales where atomic-scale forces take over. These forces 
then implement the interface breakup and reconnection. 
This is the case, for instance, of three-dimensional jets 
0- The pursuit of this idea in the specific context of 
rotating Hele-Shaw flows is the basic motivation of the 
present study. 

Recently, macroscopic hydrodynamic equations have 



also been used as a basic ingredient in the dynamics of 
nanojets ! Q, or in very thin membranes, where the 
effects of Van der Waals forces must be taken into account 
. These are examples of nonlinear processes where 
the macroscopic dynamics coupled with molecular forces 
or thermal noise lead to pinch-off phenomena. 

In the narrow gap between the two parallel glass plates 
of a Hele-Shaw cell, the scale where the effective two- 
dimensional "macroscopic" description fails is basically 
the thickness of the cell gap [7|. Although this cutoff 
is much larger (typically, of the order of the millime- 
ter) than the microscopic cutoff of the hydrodynamic de- 
scription of 3d jets, it still makes sense to investigate to 
what extent the 2d effective dynamics in a Hele-Shaw cell 
(playing the role of the "macroscopic" dynamics) lead to 
the spontaneous occurrence of finite-time singularities. 
The interface recombination will certainly introduce new 
physics (not contained in the ordinary Hele-Shaw equa- 
tions) which could easily depend on details that may not 
be universal, such as wetting conditions, 3d structure of 
the meniscus, contact line motion, etc.; the precise time 
of the pinch-off may also depend on these details. How- 
ever, the underlying idea is that the effective cutoff and 
the details of the actual 3d pinch-off can in principle be 
arbitrarily reduced, for instance changing the cell gap. 
So, ultimately, the relevant question is what the pre- 
diction of the effective 2d (Hele-Shaw) dynamics is: If 
the 2d dynamics lead themselves to finite-time pinch-off, 
this would certainly assure the existence of finite-time 
pinch-off in the real experiments, and would also define 
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an upper bound for the actual experimental pinch-off. If, 
on the contrary, no finite-time pinch-off is obtained for 
the Hele-Shaw dynamics, it will be relevant to study the 
evolution of a thin neck as predicted by the effective 2d 
dynamics, in order to elucidate, in each particular case, 
whether actual pinch-off would be expected invoking the 
additional 3d effects. In this context, the simulations 
performed within a phase-field scheme are particularly 
appropriate, because they incorporate a natural cutoff 
(the interface thickness) which controls the actual pinch- 
off and which, like the cell gap in Hele-Shaw flows, can 
be modified at will to postpone the pinch-off of a narrow 
fluid filament. 

It is well known that the effective two-dimensional 
Hele-Shaw dynamics can lead per se to pinch-off. In fact, 
it was for this type of dynamics that finite-time pinch- 
off singularities in a hydrodynamic description were first 
found [a As opposed to the 3d case of a cylin- 
drical interface, a 2d fluid filament (two parallel inter- 
faces) is stable to the Rayleigh criterion: The total in- 
terface length (at constant area) is larger for any per- 
turbation around straight interfaces. However, surface 
tension alone has been shown to drive a configuration of 
two droplets of fluid connected by a neck to finite-time 
pinch-off in two-dimensional simulations, for certain ini- 
tial conditions 0, |n| . A straight filament of fluid can 
also be made to pinch at infinite time (and sometimes 
also at finite time), with very specific boundary condi- 
tions 0, 0, 0, 0, [3 ■ However, two-dimensional pinch- 
off in unprepared situations emerging spontaneously in 
externally driven dynamics, such as in traditional vis- 
cous fingering experiments, has not been addressed. One 
of our main goals is to account for the spontaneous ap- 
proach to pinch-off singularities often observed while a 
highly nonlinear pattern develops following a morpho- 
logical instability, in particular for rotating Hele-Shaw 
flows. Pattern formation and pinch-off singularities can 
indeed be intimately related, as it will become apparent 
in the case studied here. 

Another issue that it is worth addressing is the role 
of the viscosity contrast or Atwood ratio A = (jUi n — 
Mout)/(A 1 in + Aw) (where /ij n and ^ out are the viscosities 
of the inner and outer fluids respectively) in the pinch- 
off process. It is known that this parameter has a strong 
influence on the interface shape near pinching in three 
dimensions 0, but only recently has it been remarked 
that the most studied limit A — ► 1 might be a very spe- 
cial case also in three dimensions [l^. I n effectively two 
dimensions, experiments performed in this limit, with air 
displacing a liquid in a channel geometry, show that the 
fingers formed due to the morphological instability com- 
pete until a single finger is left 0] ; the neck of the tran- 
sient and final finger(s) does not even pinch. In con- 
trast, when a denser liquid is put on top of a less dense 
one of similar viscosity (typically A ~0-0.5) in the same 
but tilted channel, fingers do not compete, but elongate 
to form thin filaments with a droplet at their tip |17| . 
which can indeed pinch off [18j. One possible scenario 



to explain this different tendency to pinch-off might con- 
sist in relating the absence of competition for low A with 
the formation of long filaments, which could then pinch 
by mechanisms not qualitati vely different from those of 
previous studies |8L l9l flOt [llL [13L [14^ . 

We have performed experiments in a rotating Hele- 
Shaw cell, where a drop of the more dense fluid is placed 
at the center and the instability is driven by the cen- 
trifugal force. The idea is to favor the formation of long 
(radial) filaments by the use of a driving force which in- 
creases as the fingers become longer. The natural ques- 
tion to pose is then whether these long filaments will also 
pinch for high viscosity contrasts. 

Long fingers are indeed formed and they stretch and 
narrow for any viscosity contrast A. However, for 
high values of A we only observe pinch-off sporadically, 
whereas we find it systematically for low ones. This rules 
out the above scenario. It rather suggests that the pre- 
cise nonlinear dynamics close to pinch-off, and not only 
the overall morphology, is also dependent on the viscosity 
contrast A. This dependence could either be contained in 
the (2d) Hele-Shaw dynamics or, instead, enter through 
different wetting properties for air-oil (high A) and oil- 
oil (low A) interfaces when the width of the filaments 
becomes comparable to the cell gap in our experiments. 

We have run phase-field simulations of the two- 
dimensional dynamics which lead to pinch-off singular- 
ities for low, but not high viscosity contrast. Although 
we can only track the first stages of the approach, the ma- 
jor dependence of the pinching dynamics on the viscosity 
contrast thus seems to be contained in the Hele-Shaw 
equations. In particular, our simulations indicate that 
finite-time pinch-off could be contained (at least) in low- 
A Hele-Shaw dynamics for more general settings than 
those originally studied in Refs. 0M ES IH 13 111 • 
Therefore, it is plausible that the effect of the cutoff (here, 
the cell gap) , although entering the problem earlier than 
in usual three-dimensional cases, can still be regarded 
as an implementation of pinch-off and reconnection, not 
necessarily affecting significantly the approach to pinch- 
off, which is governed by the Hele-Shaw dynamics. 

The rest of the paper is organized as follows: In Sec. [H] 
we present the experimental results, with an indication of 
scaling behavior. This is then explained in the following 
section by simple theoretical arguments valid asymptot- 
ically. In Sec. II VI we test this basic picture by numerical 
simulations of the Hele-Shaw equations. These confirm 
theory and experiments to agree for high viscosity con- 
trasts, whereas puzzles remain for low ones. One of these 
can be solved by rethinking the experimental time scale 
as done in Sec. [3 while another requires a refinement 
of the theory which is performed in Sec. IVII The results 
are summarized and further discussed in Sec. IVIII 
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II. EXPERIMENTAL RESULTS 

We have performed two sets of experiments in a rotat- 
ing Hele-Shaw cell: one for high and one for low viscosity 
contrast. These two types of experiments were already 
performed for other purposes in 0] and [2(j; we repeat 
them here to address the scaling behaviors and the pres- 
ence of pinch-off as functions of viscosity contrast. Our 
cell consists of two horizontal glass plates, 7 or 10 mm 
thick and 390 mm in diameter, separated by different 
spacers of height b = 0.5,0.7,0.8, 1.0, 1.4 mm. The cell 
is rotated around a vertical axis intersecting its center, 
with a controlled frequency tt, and the interface shape is 
recorded with a CCD camera. 

In a first set of experiments, a silicone oil of viscosity 
p in = 50cp and density p- m = 975 ± 10kg/m 3 at 20°C 
and air are used as inner an outer fluids, respectively, in 
a prewetted cell |2l|. This yields A = 1 and Ap = p in — 
p out ~ 0.98g/cm 3 , and an interfacial tension a — 20.7 
mN/m. In a second set of experiments, we use another 
silicone oil (/i; n = 530 ± 50cp, p- m = 975 ± 10kg/m 3 ) as 
inner, and a vaseline oil (p ut = 190±50cp, p out = 875 ± 
10kg/m 3 ) as outer fluid. In this case A = 0.45 ± 0.05, 
Ap ~ 0.10 g/cm 3 , a = 1.8 ± 0.7mN/m. Mean values are 
given at 20° C, and the uncertainty in the viscosities (and 
their contrast A) accounts for temperature variations. 

Figures ^a,b) and ^c,d) show the typical patterns 
formed at two different stages, for high (A = 1) and low 
(A ~0.45) viscosity contrast, respectively. Although the 
nonlinear dynamics depend significantly on the viscosity 
contrast, as we can see, let us focus for the moment on 
those features which are generic to any value of A: After 
some latency time, the linear instability leads to small 
undulations on the initial circle of radius Rq . These un- 
dulations grow into radial fingers of silicone oil. As they 
are further centrifuged, the fingers develop overhangs; 
their tips evolve into differentiated droplets, whereas the 
overhang regions themselves narrow and stretch to be- 
come thin filaments, which continue to stretch and nar- 
row. In parallel with this, incoming fingers of less dense 
fluid advance towards the cell center. However, they do 
not develop thin filaments nor droplets, and they slow 
down as they approach the cell center. 

The filament thinning and the dynamics of the incom- 
ing fingers are rather difficult to measure accurately and 
not very sharply defined experimentally. In contrast, the 
stretching of an individual filament, or, more precisely, 
the radial coordinate of the droplet at its tip, R, is more 
readily accessible. The squares (A = 1) and the circles 
(A ~ 0.45) in Fig. [21 indicate this droplet position (even 
after pinch-off) in units of the radius Rq of the initial 
circle and in log scale, as a function of time t in units of 
the time scale 

_ 12(p in + p out ) 

WAp ■ ( ' 

The open squares and circles indicate the evolution of the 
furthest droplet in Figs. nja,b) and^c,d) respectively. 




FIG. 1: Pattern evolution for A = 1, b = 0.5mm, Q = 
120rev/min and Ro = 50 mm (a,b), and A ~ 0.45, SI = 
180rev/min and Ro = 38 mm (c,d). Snapshots (3Ro x 3Ro) 
16.5s (a) and 22.5s (b), and 122s (c) and 158s (d) after be- 
ginning rotation. 



As we can see, the growth tends to be roughly ex- 
ponential, at least for R/Rq > 2. (Note that R/Rq < 2 
corresponds rather to relatively small fingers, rather than 
droplets at the tip of thin filaments). It seems reasonable 
to measure a growth rate — the slope of linear fits (solid 
lines) in Fig. |2 this yields to = 0.8-1.0 with good linear- 
ity beyond R = 2Rq for A = 1, and to = 1.25 — 1.5 
for A ~ 0.45. A ~ 0.45 curves tend to bend a bit 
more, but this eventual bending is uncorrelated with the 
droplet pinch-off: To check this, we have performed an 
experiment with an isolated, off-center circular droplet; 
it roughly keeps its circular shape and runs away at a 
very similar rate (to = 1.45, triangles). 

The most striking difference between high (A = 1) 
and low (A ~ 0.45) viscosity contrast is the system- 
atic droplet pinch-off observed for low, but not high A. 
For the latter, filaments keep on growing and stretching, 
while reaching a width comparable to the cell gap, that 
is close to the natural cutoff. While we have not spanned 
the whole range of A, all previous evidence shows that 
the so-called low contrast behavior is qualitatively un- 
changed for most of the range of A, and only very close 
to A = 1 significant differences are found. This has been 
recently discussed in Refs. flii, |22|| . 

Another difference seems to be the growth rate of 
droplets (connected or not to filaments). The uncertain- 
ties in the viscosities are not sufficient to account for it: 
droplets seem to be centrifuged faster in dimensionless 
time for low than for high viscosity contrasts. Our dis- 
cussion on this is postponed to Sec. |V] 
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FIG. 2: Droplet kinetics: Evolution (as a function of scaled 
time) of the radial coordinate of (i) the center of mass of 
the drop at the tip of three different filaments for A ~ 0.45, 
before and after pinch-off (circles); (ii) the center of mass of 
an isolated drop for A ~ 0.45 (triangles); (iii) the tip of three 
different filaments for A = I (squares). Circles and triangle 
have been translated to earlier times for comparison. Open 
symbols correspond to the most excentric drops in Fig. m 
are the slopes of linear fits. 



III. THEORETICAL SCENARIO 

We try to link here by continuity arguments the fila- 
ment thinning and stretching and the evolution of incom- 
ing fingers, and calculate the various growth and decay 
rates in a certain approximation. Because the dynam- 
ics of droplets at the tip of filaments should affect their 
stretching, we also construct an exact time-dependent 
solution for an off-center isolated droplet. Finally, we ex- 
plore the existence of stationary-shape solutions as possi- 
ble asymptotic attractors suggested by the slowing down 
of incoming fingers. 



A. Scalings 

Fluid incomprcssibility imposes fluid area conserva- 
tion, which intuitively forces filaments to stretch and/or 
feed the droplet at their tip as they narrow, and fingers 
of less dense fluid to grow inwards as mass is transported 
outwards through the filaments. Area conservation hence 
seems an important feature of the dynamics. In partic- 
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FIG. 3: Sketch of a filament of fluid with a droplet at its tip. 



ular, "star-fish" patterns have been observed in an ad 
hoc geometrical model for the normal velocity of a front 
with local mass conservation on the interface (+ a con- 
stant flux) 23] . However, incompressibility also generally 
leads to non-locality, which will turn out to be crucial for 
low viscosity contrasts and for the presence of pinch-off 
singularities. 

Let us first address the filament dynamics by means of 
a continuity equation in the overall direction of the fila- 
ment x, which reflects the incompressibility of the inner 
fluid: 

d t h + d x (hv x ) = 0, (2) 

where h[x, t) is the local y coordinate of either the upper 
or the lower interfaces, and v x is the x component of the 
fluid velocity inside the filament (see Fig. 0J). 

Motivated by the experimental observations, we look 
for a scaling solution of the form 

h(x,t) = L(t)h [xL{t)], (3) 

where L(t) is a dimensionless scaling factor. Substituting 
this ansatz into Eq. we find 

d t L h d x v x + v x d x h 

T~ = 1 i Q~T = K > ( 4 ) 

L no + xa x no 

where K is independent of x. 

If v x happens to be either (i) a constant along the 
filament or (ii) proportional to x, v x — Kx, a per- 
fectly straight filament (ho independent of x) satisfies 
Eq. (@J and hence stays straight. Moreover, if and only 
if v x = Kx, any function ho(x) satisfies Eq. J3J and 
preserves its form. Interestingly, we will see that, for a 
straight filament, (i) corresponds to the standard chan- 
nel geometry, and that our setup mostly falls in case (ii) . 
In contrast, the radial case with injection, for example, 
cannot exhibit this type of self-similar solutions. 

We find that, whenever the above scaling ansatz is 
valid, and further assuming K to be constant in time, 
then L = e~ Kt . The initial condition ho(x) thus narrows 
exponentially 

h = h e- Kt (5) 

and stretches 

R n = R n {t = 0)e Kt (6) 
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(where R n = x | dtho(x,t) = 0) with a time constant 
set by the inner fluid velocity. Figure schematically 
shows the in- and outgoing flow in this case. Eq. © 
is a statement that a reference point on the interface is 
advected with the filament as this is stretched. In partic- 
ular, if applied to the neck of the filament just before the 
droplet as represented in Fig. |3 it would predict the to- 
tal filament length to grow with the same time constant 
K with which its width decays. However, note that the 
outgoing flux of inner fluid at the neck can go into the 
droplet or contribute to the filament stretching. If some 
fraction of it actually goes into the droplet (as observed 
in the experiments), the scaling of the radial coordinate 
of the neck will be delayed with respect to the thinning 
rate K. 

Finally, in order to relate this filament scaling with 
the incoming fingers, consider a reference circle of radius 
i? re f intersecting all filaments near their root, where they 
begin to look straight. The outwards flux of denser fluid 
through the circle and inside each filament is 2h(t)v x (x — 
R Ie f), and assuming again that v x (R le f) does not depend 
on time, we find it to scale as e~ Kt . In general, the area 
of denser fluid enclosed by this reference circle a should 
hence asymptotically decay with the same rate, 

a = a(t -> oo) + [a(t = 0) - a(t -> oo)] e~ Kt . (7) 

A similar slowing down can then be expected for the 
radial position of the tips of identical incoming fingers. 



B. Exponential behavior and time constant 

We now show that, for a straight filament, v x is linear 
in x as we anticipated, and we find the time constant K, 
which we have seen to set the time scale of the various 
exponential behaviors discussed so far. For that purpose, 
let us compute the tangential velocity jump across the 
interface. 

The 2d velocity field Ui is proportional to the gradients 
of the pressure field pi (Darcy's law), where the subscript 
i labels the inner (i =in) and outer (i =out) fluids. For 
a rotating cell 0] 

^ = ~^Pt-Pt^ri), (8) 

where r is distance to the rotation axis and f is a unit 
radial vector pointing outwards. Using the standard 
boundary condition of Laplace law for the capillary pres- 
sure jump across the interface, p\- ln — p\ ou t = ck, where 
the bar | stands for the limit value when approaching the 
interface from either side and k is the interface curvature 
in the cell plane, we obtain 

b 2 

s ■ (u|out - u|in) = ; r (o-d s K - Apf2 2 rf • s) 

6\J4jx + //out) 

+A§- (u| in +u|out), (9) 



where s is a coordinate tangential to the interface and s 
the unit vector in that direction. 

In our sketch (Fig. EJ), we identify v x ~ u\- m ■ s, which 
is a function of the local interface geometry and (1 — 
^)u| ut • s, the solution of a non-local problem. To keep 
the problem local, we now assume that the outer fluid 
approaches the filament normal to it, w| ut -L s (Note 
that this assumption is not necessary for A = 1). 

For a straight (k = 0) radial (x \\ f) filament, we obtain 
the anticipated v x — Kx and the time constant 

K _ b 2 Aptf 2 1 = b 2 A P tf 

Mout 

)(1 + A) 1 + At* 12/i in v ; 

[case (ii) above]. One might wonder how a scale- free so- 
lution is possible without requiring the filament width to 
be much smaller than any other length scale, as it is le- 
gitimate to do close enough to pinch-off pj. The answer 
is that the assumption of flatness has removed the capil- 
lary length (proportional to a) from the local problem; it 
still sets, together with the radius of the initial condition 
Ro, the size of e.g. the droplet at the tip of the filaments, 
but no external length scale can be felt when dropping 
the nonlocal term (1 — A)w| out • s. This allows for scal- 
ing solutions to exist well before pinch-off for A = 1. It 
is this approach to pinch-off, rather than the pinch-off 
phenomenon itself, which will be further studied in this 
article. 

For a tangential filament, K = and the filament 
is stationary [case (i) above], which is consistent with 
our observations for the tangential segments of some fil- 
aments. These segments are found where a tip split (for 
A = 0, see Fig. ^c,d)) or at the point where an incom- 
ing finger was overcome and stopped by its neighbors (for 
A=l, see Fig. Eta.b)). 

For the case of viscous fingering in a channel, the linear 
increase in the centrifugal force is replaced by a constant 
gravity (or equivalent injection). We can easily track this 
difference to yield a fluid velocity inside a filament par- 
allel to the gravity or injection direction v± independent 
of x. The filament width will again be stationary [case 
(i) above]. We see by comparison that rotation not only 
ensures the formation of long filaments as we expected, 
but also brings the two interfaces close together. 

For not completely straight filaments h(x), v x ^ Kx. 
Therefore, strictly speaking, the scaling solution Eq. (J3J) 
does not hold for any ho(x), but only for a straight fil- 
ament. It remains an approximate solution for slightly 
curved filaments (see Sec. IVI|) . but it makes little sense 
for the droplets at the tip of the filaments. This raises the 
issue of dynamics of the droplet, which could potentially 
affect the filament stretching. 

C. Isolated droplet 

The velocity found for the filament stretching has an 
elementary interpretation: In the overdamped limit as- 
sociated to Darcy's law, all forces balance. The inner 
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fluid velocity v x will hence be that which makes friction 
balance the other external forces, namely the centrifugal 
force and the force exerted by the outer fluid through its 
pressure gradients. The latter can be derived from an- 
other force balance, now in the outer fluid. There, the 
viscous forces vanish in the x direction with our assump- 
tion u out -L s ~ x, or even without this assumption for 
A = 1, since it corresponds to p ou t — 0. Thanks to this, 
the other forces in the outer fluid (centrifugal and pres- 
sure gradients) balance in this direction. Thus, back in 
the inner fluid, friction, \2p, m v x /b 2 , only needs to com- 
pensate for centrifugal forces, ApVl 2 x. This gives rise to 
the velocity found, v x — Kx with K given by Eq. (|10|l . 

For a droplet (a closed interface), in contrast, the outer 
fluid needs to be moved in the x direction. The outer 
viscosity and velocity field enter the force balance, and 
v x < Kx, where v x is now the velocity of the center 
of mass of the droplet, and the equality holds only for 
A = 1. Except for A = 1, it is clear that the filament 
stretching will indeed be affected by this lower droplet 
velocity. Furthermore, v x will depend on the outer flow 
field, and hence on the droplet shape. 

For a slightly off-center circular droplet, v x is known: 
Such a droplet corresponds to the mode 1 in a linear 
stability analysis of a circular, centered drop. Its growth 
rate is K' = 1/t* = b 2 Apfl 2 /[l2(p in + p out )] [TU, where 
we readily see that it is the sum of viscosities, and not 
only the inner one, what counts. 

We now compute the velocity for a circular droplet 
d s n = rigorously. Preserving the circular shape re- 
quires that Uin = R x everywhere (R is distance from 
the rotation axis to the center of the droplet), so that 
the circle be just translated as a whole. We express the 
unit vectors tangential (s) and normal (h) to the inter- 
face and the radial vector r connecting the rotation axis 
and a point on the interface in polar coordinates (/?, a) 
with respect to the center of the droplet (see Fig. |3J), 
and compute s ■ u\- m = —Rsma and r • s = — Rsma. 
For A = 1, Eq. reduces to a simple geometric re- 
lation between the two. Interestingly, like the filament 
scaling, this is verified either for a channel under injec- 
tion or gravity (R is then a constant) or for a rotating 
cell as presented here, R = K'R, but not for a general 
driving force or geometry. We hence find that a circular 
off-center droplet is a time-dependent solution in a chan- 
nel or rotating cell for A = 1 as was previously found 
using conformal mapping techniques in p4j |. 

For A < 1, the outer flow enters the problem. If a cir- 
cular droplet is to remain a solution, it has to match the 
linear regime. At that stage, we will still have R = K'R, 
regardless of the viscosity contrast. Substituting this into 
Eq. ©, we find s • u| ut = — s ■ u- m = i?sina. Con- 
tinuity of normal velocities across the interface gives us 
n ■ "lout = i?cosa, and we thus find the outer velocity 
on the interface to be w| ut = R(x cos 2a + y sin 2a). We 
then propose the outer velocity field to be the product 
of its value on the interface and a function of p only. We 
find that this can indeed fulfill incompressibility if the 



droplet is isolated, i.e., if there are no other boundary 
conditions to be satisfied elsewhere. For the channel ge- 
ometry, this corresponds to the limit in which the droplet 
is much smaller than the distance between the walls. The 
complete solution reads 

Win = Rx, (H a ) 
Mout = R(po/p) 2 (x cos 2a + j/ sin 2a), (lib) 
R = R(t = G)e K>t . (11c) 

The corresponding streamlines in the outer fluid p/po oc 
sin(a), are circles of different radii tangential to y = 
and passing through the center of the droplet, although 
they change into horizontal lines inside it (see Fig. 
bottom right). 

Hence, an off-center circle turns out to be an exact 
time-dependent solution for any viscosity contrast and 
independent of it, even beyond the linear regime. For 
A = 1, K' = K and the filament stretching and the 
droplet runaway coincide regardless of the droplet shape 
and whether it is near a filament or not. This single time 
constant should hence be robust, and indeed an only 
slightly lower value [(0.8-1.0)/i*] is systematically ob- 
served in our experiments (squares in Fig. |2J). For A < 1, 
in contrast, K' < K, and an isolated, circular droplet has 
a time constant down to twice smaller (for A = 0) than 
that of the filament stretching. For a droplet connected 
to the tip of a filament (neither isolated nor circular), we 
do not know its dynamics rigorously. Our experiments 
are puzzling in this respect: Isolated droplets (triangles 
in Fig. |2l and those connected to a filament (circles) dis- 
play the same growth rate; consistency with the theory 
then requires it to be K' (the same than for A = 1), since 
the isolated droplet solution is exact. However, we seem 
to find a growth rate of K [(1.25-1. 5)/t*], and, in any 
case, higher than the measured value for A = 1. Our 
simulations will help clarify this point. 

D. Stationary-shape solution 

The various scaling behaviors derived so far point at an 
asymptotic pattern enclosing a finite area [Eq. Q; see 
Fig. 0] with radial, semi-infinite [Eq. ©] and straight 
[Eqs. © and ©] filaments of zero width [Eq. ©] (see 
Fig. 0} with a droplet at their tip. 

In a stationary state, Ui n = u ou t = by definition, and 
the interface shape must balance exactly surface tension 
with the centrifugal force in Eq. ©: 

ad s K = Aptfr ■ s. (12) 

Noting that s — dr/ds, where s is arclength along the 
interface, we can integrate the above equation once to 
find 

k(s) = k{s = 0) + A P n 2 /{2a) [r 2 (s) - r 2 (s = 0)]. (13) 

We see that the curvature of the interface k increases 
monotonically with the distance r to the rotation axis. 
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FIG. 4: Stationary solutions in a rotating H-S cell with four- 
and six-fold symmetry. 



In general, the interface will go from convex to concave 
as r is increased. The solutions of Eq. (|12ll which do 
change concavity do not resemble our observed patterns, 
but have a petal shape and are unstable [2{| . We hence 
turn to purely convex stationary solutions of Eq. (|12fl . 
up to the point where k = 0. These will represent the 
asymptotic incoming fingers. To connect them with the 
radial, straight filaments, we impose locally radial and 
straight ends at an arbitrary distance r(s = 0) from the 
rotation axis as boundary conditions. Any other connec- 
tion would cause d s n in Eq. ljl"2*|) to diverge. 

For given experimental parameters Apfl 2 /a, this com- 
pletely determines the solution, an arch covering a certain 
angle. This angle is a continuous function of the only di- 
mensionless ratio we can construct, ApQ 2 r(s = 0) 3 /er, 
and hence only for discrete values of this ratio is the an- 
gle a submultiple of 2tt. This then allows one to connect 
an integer number of equal arches to form an apparently 
closed interface, as we show in Fig. 0] This could not 
be a stationary solution on its own, since actually clos- 
ing the interface would imply that the arches meet in 
cusps. We continue each arch end by a straight, semi- 
infinite radial line, replacing each cusp with an infinite, 
straight filament of zero width. These filaments present 
no flux inside (finite velocity but zero width) nor induce 
any flow outside. They are hence compatible with sta- 
tionary arches. Finally, the interface is closed by placing 
a droplet of arbitrary size at the tip of each infinite fila- 
ment. Although a moving droplet does induce an outside 
flow around it, this finite disturbance vanishes near the 
center for an infinite filament. 

Thus, we have constructed a solution which shape is 
stationary at finite distances. Its physical relevance, how- 
ever, remains to be checked. The slowing down of in- 
coming fingers derived above and the flattening of an 
initially curved filament h{x) [Eqs. © and ©] are sug- 
gestive that such a solution can indeed be an attractor of 
the dynamics for initial conditions with the appropriate 
n-fold symmetry, and the numerical simulations in next 
section and the results of Sec. IVII further support this 
idea. 

A related issue is the selection of the length scale 
r(s = 0) of the stationary arches. For fixed boundary 
conditions, only discrete values of the dimensionless ra- 
tio Apil 2 r(s = 0) 3 /<7 allow equal arches to meet without 



crossing. Therefore, for given experimental parameters, 
the stationary-shape solutions only exist for discrete val- 
ues of r(s = 0) (note that unequal arches cannot meet at 
a single distance r(s = 0) from the center). This sets up 
a selection problem for the values of r(s = 0). 

Dynamically, such a solution will be approached with 
a certain number of fingers. This number may not be 
unambiguously defined for unequal finger configurations, 
and it is indeed a dynamic variable for high viscosity con- 
trasts. Furthermore, the initial condition is generically 
nonsymmetric, with some fingers competing and advanc- 
ing over others. As a consequence, real experiments and 
numerical simulations can evolve towards configurations 
close to the above stationary-shape solutions during a 
transient but will eventually depart from them. How- 
ever, for initially equal fingers as those one can produce 
in simulations (see next section), the dynamics conserve 
the number of fingers regardless of the viscosity contrast. 
Hence, only the stationary-shape solution of length scale 
r(s = 0) with the number of fingers of the initial condi- 
tion can be approached. This then unambiguously selects 
the scale r(s = 0) and the area enclosed by the stationary 
interface a(t — > oo) (Eq. 0). 

We thus see that the combination of the scaling solu- 
tions and the stationary state derived could specify the 
asymptotic dynamics of equal fingers and the ultimate in- 
terface shape. As for its relevance for unequal fingers, we 
point out two facts: First, we find reasonable agreement 
between the exponential rate predicted by this simple 
theory and that measured experimentally for the radial 
position of the droplets at the tip of the filaments. Sec- 
ond, the existence of the above symmetrical stationary 
solutions could explain the experimental observation of 
plateaux in the evolution of the area of inner fluid con- 
tained in a reference circle, as a function of time in Ref. 
|Tg| . In fact, in those experiments the interface could suc- 
cessively visit the neighborhood of stationary states for 
different numbers of fingers, and correspondingly slow 
down the dynamics close to them until some event (typ- 
ically a pinch-off) would trigger again the dynamics. 



IV. NUMERICAL TESTS 

We test the basic picture presented above by phase- 
field simulations of the full Hele-Shaw dynamics. They 
are compared to experimental results and to a numerical 
solution for the stationary arches. 



A. Method 

First we adimensionalize the problem: lengths will be 
measured in units of the initial drop radius R$, and time 
in units of t* as given by Eq. QJ. All quantities ap- 
pearing below will be dimensionless. The scaling laws 
in the previous section can be casted into dimensionless 
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form just by considering all quantities appearing there 
dimensionless and redefining K = 2/(1 + A) and K' = 1. 

The standard Hele-Shaw dynamics is defined by the in- 
compressibility condition already used before (V-tt, = 0), 
together with the continuity of normal velocities across 
the interface and the tangential velocity jump [Eq. JJjJ] 
as boundary conditions. Thanks to incompressibility, we 
can define a stream function u — V x ipz (where z is the 
direction perpendicular to the plates), which is contin- 
uous across the interface. We thus obtain the following 
governing equations [2^, : 



(14a) 
(14b) 



V 2 ip = o, 

d s ^\in = d s 1p\ out = -V n 



dn l 



On lm ' 



(14c) 



where the first expresses that the flow is irrotational in 
each fluid, the second is the continuity of normal veloci- 
ties, and the third, their tangential jump. The magnitude 
of this jump is the strength T of a divergent vorticity 
peaked on the interface, T = 7 + A(d n ip\i n + d n %jj\ ou t), 
where j/2 = (BV n — rf) ■ s is its local part. We see that 
the dynamics are controlled only by two dimensionless 
parameters: the viscosity contrast A and the dimension- 
less surface tension B = a j [Apfl 2 R^] (the ratio of stabi- 
lizing to destabilizing forces). Remarkably, none of the 
scalings derived in the previous section depends on B. 

Simulations were run using a phase field model pre- 
sented in [2j| and extensively tested in [2^| : 



1 1 



e^- = V> + AV • (0V</>) + 7=7(^(1 - ),(15a) 



dt 



,89 



e-^ = f(9) + e 2 V 



Of 



e 2V2 



l 9 + e 2 K (9)\V9\ 



+ e 2 z ■ (W x V0), 



(15b) 



where 9 is the phase field, an auxiliary field distinguishing 
between the two fluids, f(9) = 9(1 - 9 2 ), ^ = s(9) ■ 

(BWk(9) + y) and k(6) = -V • h(9), with n(9) = 4^ 

and s(8) = h(9) x z. 

Apart from the physical control parameters A and B, 
the dynamics in this model also depend on an artificial 
interface thickness e and a relaxation time for the stream 
function (which is diffusive, not Laplacian) e. In the limit 
e, e — ► 0, the dynamics are strictly those of Eqs. (|14al - 
I14c|) . For finite values of e,e, an error bound is guaran- 
teed for any given magnitude by conveniently decreasing 
e,e|23. 

A very interesting feature of this type of model is that 
both fluids and the interface between them are treated 
as bulk. One consequence of this is that the model is 
well behaved as two interfaces break up and reconnect, 
which enables us to study the dynamics after pinch-off 
in a very natural way. This has been demonstrated very 



recently |30| for the model used here |28|, |29| reformu- 
lated in terms of the velocity vector. Let us add a note 
of caution: When two interfaces approach to distances 
comparable to their thicknesses e (more precisely, below 
5e from numerical tests) they attract each other. This 
triggers the pinch-off. One could think of this effect as a 
phenomenological "pinch rule" . We nevertheless do not 
pretend the phase-field dynamics near pinch-off to rep- 
resent the effect of the third dimension nor the eventual 
breakdown of a hydrodynamic description. Our approach 
is to simulate only the bidimcnsional Hele-Shaw dynam- 
ics down to distances of 0(e), and continue to simulate 
them once the reconnected interfaces have separated a 
distance of the same order, very much as in the experi- 
ments, where the Hele-Shaw dynamics can be considered 
to be accurate before and after the two interfaces were 
separated a distance of 0(b). The crossover from the 
Hele-Shaw dynamics to interface interaction via their fi- 
nite thicknesses is very abrupt, enabling us to clearly 
separate one from the other. 

For reliable quantitative comparison with theory and 
experiments, we explicitly check convergence of the time 
constant K in e, e, since e conveys some finite diffusion 
time to the flow, and e delays the interface advance with 
respect to the normal fluid velocities [28|. We thus es- 
tablish exponential rate values for the e, e — > limit, i.e., 
for the Hele-Shaw dynamics, and bounds for e, e to ob- 
tain those values in practice. For e, a value of e = 0.5 
turns out to suffice, and it is used in the following unless 
otherwise stated. Convergence in e is discussed case by 
case. 

We numerically integrate the above model with an Eu- 
ler scheme and centered differences. The time step dt is 
taken to be close to the stability limit; dt = 0.2edx 2 un- 
less otherwise stated, where dx is the mesh size. Conver- 
gence of the solution is also tested changing dx (dx = e 
generally suffices as shown in |2{|). Our initial condition 
is some perturbation of a centered circle of unit radius, 
i.e., r = 1 + Ar(ip) where r, ip are polar coordinates 
with respect to the rotation axis and Ar(ip) < 1. In the 
next two subsections, we consider only identical fingers 
(n-fold symmetry), Ar(<p) = (2n /n)qcos(nip), where q is 
the amplitude to wavelength ratio of the perturbation, 
and n is the (integer) number of fingers. Unless other- 
wise stated, we use a dimensionless surface tension B for 
which n is the most unstable mode in the linear regime 
and q = 0.05. 



B. Filament thinning 

We test here the filament thinning scaling by looking 
at a dumbbell-shaped pattern: two droplets connected 
by a filament (n = 2, see Fig. This case has two 
advantages: a width can be unambiguously defined as 
that at the axis of rotation (the midpoint between the two 
droplets); on the other hand, there are no neighboring 
filaments to influence its dynamics. 
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FIG. 5: Filament thinning: Initial condition (top, q = 0.05) 
and one later interface for B = 0.09, e — 0.008, dx = e/2 and 
A = 1 (middle) or A = (bottom). 



Figure displays the patterns thus obtained at the 
end of our run (A = 1) or just after the first pinch-off 
(A = 0). In the last case, new pinch-off events quickly 
further shorten the central filament from its ends. We 
will refer to this phenomenon as "pearling". The two 
upper curves in Fig. [fj] (B = 0.09) show the evolution 
of the filament width at the midpoint for the whole run 
(A = 1) or up to the relaxation of the last central segment 
in a round shape, which makes the width increase again 
(A = 0). The latter occurs when pearling reaches the 
rotation axis, shortly after the first pinch-off, indicated 
by the vertical lines. We see that the asymptotic decay 
is grossly exponential for A — 1 as predicted, although 
it is only observed for roughly a decade. For A = 0, the 
filament thinning is strictly not exponential since the ef- 
fective exponential rate varies slowly in time. It should 
be noticed that it is the flatness, not the width of the 
filament, what determines the goodness of the approx- 
imation of last section. One might think that pinch-off 
prevents the formation of long and hence straight enough 
filaments for A = 0, but, indeed, just before it the A = 
filament is just as straight as the A = 1 one at the same 
time; the A = 1 filament already exhibits a clear expo- 
nential behavior, whereas the A = one does not. 

A closer look reveals that the thinning rate is well re- 
laxed to its asymptotic value for the last third of the run 
for A — 1, with a value of K = 0.98 consistent with the 
predicted K, but not for A = because pearling reaches 
the center when the rate begins to approach (from above) 
a value of K = 2. No disturbance is observed from the 
previous pinch-off events. Runs for A = 0.5 and A = 0.8 
(not shown) are also in reasonable agreement with the 
law K = 2/(1 + A). 

One might blame the high dimensionless surface ten- 
sion B — 0.09 used or the finite interface thickness e or 
relaxation time e in the model for either the slow relax- 
ation of the exponential rate or the pinch-off itself in the 
low-^4 case, since B enters the problem as long as the fil- 
ament is not perfectly straight. The lower curves in Fig. 
[5] for a lower dimensionless surface tension, B = 0.01, 
and at different values of e and dx/e rule out both possi- 
bilities. Using a more deformed circle (q = 0.2), to avoid 
the appearance of other modes which arc linearly unsta- 
ble, we can confirm the results obtained for B = 0.09. 
For high viscosity contrast, the good decay even for the 
coarsest mesh and interface profile allows to establish the 




1 pinch-off 2 pinch-off 3 

t/t* 



FIG. 6: Filament thinning: Width at rotation axis (in log 
scale) vs. time. Dashed, solid and dotted line(s): dx = e, 
dx — e/2 and dx = e/4, respectively. Upper curves (B = 
0.09): e = 0.008 (runs in Fig. q = 0.05). Lower curves 
(B = 0.01, q = 0.2): three curves at e = 0.02, 0.01 and 
0.005 for dx = e and three at t = 0.02, 0.01, and 0.0067 for 
dx = e/2 shown for each value of A. For A — 0, lower curves 
correspond to lower values of e. The only dx = e/4 curve is 
for A — 1, e — 0.02. The dashed (solid) vertical lines indicate 
the interval during which the droplet at the tip pinches off for 
dx = e (e/2). 



value of the rate K unambiguously at K — 0.99— 1.00 by 
studying the convergence in e and dx. For low viscosity 
contrast runs (A = 0), the values of the rates again seem 
to relax towards the end of the runs to a value compat- 
ible with the theoretical value K = 2 although they are 
again badly defined. 

Pinch-off is systematically observed for low (A = 0, 
0.5), but not high (A = 0.8, 1) viscosity contrasts. For 
the latter, a numerical instability due to noise in the (rel- 
atively) inviscid fluid often sets in abruptly at late times, 
especially for the coarser discretization dx = e. Because 
it is so abrupt, it is easy to separate from valid data. For 
the former, the pinch-off itself cannot be considered a 
spurious effect due to a finite interface thickness, since the 
pinch-off time increases only marginally when decreasing 
the interface thickness e, and in any case stays bounded 
by the vertical lines in Fig. [S] Only for the coarsest 
interface profile can we detect that the center of the fila- 
ment anomalously "anticipates" the first pinch-off at its 
tip; for lower values of e it actually continues to narrow 
unaffected after the first pinch-off and until the subse- 
quent pearling reaches it. Interestingly, when decreasing 
the mesh size dx from dashed to solid lines, pinch-off is 
appreciably accelerated, so that it must also remain in 
the continuum limit. 

We conclude that the theoretical exponential thinning 
and the values of the exponential rates derived in last 
section are a good approximation to the observed behav- 
ior, especially for high viscosity contrasts. Deviations do 
exist and are intrinsically much stronger for low viscos- 
ity contrasts, but, because they are not less apparent for 
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FIG. 7: Approach to the stationary shape for B — 0.0213 
and e = 0.01. Left: initial (dotted) and latest (solid line) 
interfaces for A = 1, dx = 0.5e and e = 0.5. Right: Blow-up 
of central region with intermediate interfaces (solid lines) at 
constant time intervals. Thicker line: receding interface (after 
pearling) for A = 0, dx = e and 1=1. Dashed line: numerical 
solution of Eq. 11131 for the predicted stationary shape. 



a surface tension almost one order of magnitude lower, 
they cannot be due to it. 



C. Stationary-shape solution 

We now increase the number of fingers arising from the 
linear regime to n = 4, so that the incoming fingers en- 
close a finite area. This enables us to test the approach to 
the stationary-shape composite solution derived in Sec. 

inna 

Figured (left) shows the initial condition (dotted line) 
and the last interface (solid line) of our most refined 
(e = 0.01, dx — 0.5e and e = 0.5) run for A = 1, just 
before spurious pinch-off occurs. The latter is unavoid- 
able because the filament width reaches 0.05-0.07, at the 
edge of diffuse-interface attraction (~ 5e). A — runs 
pinch much faster, presumably due to a physical phe- 
nomenon and not to the finite interface thickness, since 
here, as well as for n = 2, this is not cured by decreas- 
ing e. Due to this different behavior, A = runs do 
not lead to the star-like pattern shown here, but, after a 
first pinch-off event at moderate radii, they quickly emit 
several droplets while the central drop recedes. 

The blowup in Fig. d (right) shows the last interfaces 
for A = 1 at constant time intervals (thinner solid lines) 
compared to a numerical solution of Eq. H13|) for the com- 
posite stationary pattern derived in Sec. 1111 Dl (dashed 
line). The latter seems a plausible asymptotic solution in 
the sense that it has not been overcome by the incoming 
fingers of our run, although visually these seem to asymp- 
totically approach less curved arches. For a quantitative 
estimate, we fit a slowing down of the type of Eq. (JJJ 
[i.e., we fit a(t — > oo), a(t = 0) — ait — > oo) and K] both 
to the time evolution of the area enclosed by a reference 
circle of various radii R Ie {, and to that of the radial coor- 
dinate of the incoming finger tips. We find K = 1.06 for 
the area enclosed by a circle of radius i? rc f = 1.25, where 



the filament begins to look straight, and K — 0.9 — 1 
for the incoming finger tip. Both are in good agreement 
with the theoretical value K = 1. The area enclosed by 
smaller circles decays with smaller exponents (0.93 for a 
circle of radius R rc f = 1 and 0.81 for i? re f = 0.75), but 
a worse agreement is understandable, since the filaments 
are less straight at the latest times we can explore. 

Regarding the tip of the incoming fingers and the to- 
tal enclosed area in the center, the numerical solution 
of Eq. 1|13|) predicts r m i n = 0.133 and ait — > oo) = 
0.0937 respectively. The numerical simulations here pre- 
sented show that the filament still contains a significant 
amount of inner fluid, yielding an asymptotic enclosed 
area which decreases with the radius of the reference cir- 
cle (a(t -> oo) = 0.36 for R Tc{ = 1.25, a(t -> oo) = 0.27 
for R tBi = 1 and a(t -> oo) = 0.20 for i? rcf = 0.75). 
As for the asymptotic tip position of the incoming fin- 
gers, we find r m i n = 0.2, with little sensitivity to the 
time window used. We thus seem to find that the tips 
and the area enclosed reach roughly half of those pre- 
dicted by Eq. (|13|l . This mismatch is most probably due 
to the fact that we do not yet probe close enough to the 
stationary-shape solution. This is indeed the case regard- 
ing the filaments since, as explained above, they have not 
reached its zero- width asymptotic state because the final 
area is not independent of the reference circle radius (as 
it should be if i? rc f < r(s — 0)). 

However, the evolution of the interface for A = af- 
ter pinch-off does hint that this exact solutions have an 
important effect on the dynamics. After the pearling 
process of the filaments for A = 0, the incoming fingers 
continue to approach the rotation axis and the droplets 
left at the tip of the filaments recede. When they reach 
the central region roughly delimited by the ends of the 
arches of our predicted stationary-shape solution, surface 
tension finally stops the advance of the incoming fingers. 
The whole drop left at the center relaxes towards a cir- 
cular shape. The thicker solid line in Fig. d (right) cor- 
responds to such a central shape at the time when the 
incoming fingers stop before the relaxation to a circular 
shape. Although the shape displayed is not an asymp- 
totic one, the remarkable resemblance of the tip of the 
incoming fingers and the steady-state asymptotic arches 
(dashed lines) and the fact that incoming fingers precisely 
stop (their velocity crosses zero) when approaching this 
shape is suggestive that our steady-shape solution (or 
parts of the solutions found in (23 very similar to it) 
is indeed the attractor for sets of initial conditions with 
n-fold symmetry. 



D. Unequal finger dynamics 

Finally, we address the more realistic case of unequal 
fingers. The goal here is no longer to test our theoret- 
ical results, but to directly compare with experimental 
patterns. We thus check whether the differences between 
low and high viscosity contrasts A observed in experi- 
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merits, including systematic or not pinch-off events, can 
be explained by the Hele-Shaw equations. 

In order to compare with the patterns shown in Fig. 
^ which have 15-16 filaments for A = 1 [Figs. IHa,b)] 
and 20-21 for low A ~ 0.45 [Figs. Qfod)], we use a di- 
mensionless surface tension of B = 1.03 x 10~ 3 , which 
corresponds to a most unstable mode n — 18. The idea 
is to use exactly the same initial condition and physi- 
cal and computational parameters for both low and high 
viscosity contrast A, except for A itself. Thus, any dif- 
ference we observe can only be due to the different values 
of A. 

We adopt the less refined choices dx — e and e = 1 
with an interface thickness down to e = 0.005, in order 
to avoid possible spurious pinch-off events. As already 
observed for the dumbbell-shaped pattern in Sec. IIVBI 
we find that, whenever pinch-off occurs, a lower rate of 
dx/e does not prevent it, but even accelerates it. 

As for the initial condition, if we start with a pertur- 
bation where sine and cosine modes have uniformly or 
Gaussian distributed random amplitudes up to a certain 
cut-off (here n = 40), we cannot reproduce the strong 
asymmetry of the experimental patterns for high viscos- 
ity contrast A = 1 [Figs. IHa,b)]. We traced this back 
to the fact that, in the experiments, the deviations from 
a circle first become visible in particular spots, and not 
uniformly on the interface. For the experiment of Figs. 
EIa,b), this happens roughly on three spots. We mimic 
this by multiplying our random initial condition by an ex- 
ponential envelope peaked at three uniformly distributed 
random angles, keeping Ar(ip) < 0.004. We do not at- 
tempt to truly copy the experimental initial condition, 
but to use a somehow "statistically" similar one. 

If some details look alike, these will then indicate some 
general common trends representative of the dynamics, 
as opposed to a particular initial condition. 

This initial asymmetry in the deviations from a circle 
is largely preserved throughout the whole evolution for 
high [A = 1; Figs. E[a,b)], but not low viscosity con- 
trast [A = 0, 0.5 runs; Figs. Efod) and Figs. E£e,f), 
respectively]. We hence conclude that both an asym- 
metric initial condition and a high viscosity contrast are 
necessary to obtain asymmetric patterns. This effect of 
the viscosity contrast is due to the incoming finger com- 
petition for high A, although it is less apparent here than 
in the channel geometry, due to the geometric decrease 
of available room fingers experience. 

Figure [S] is indeed the simulation analogue of Fig. ^ 
Experimental patterns in Fig. ^ were rotated at will to 
make the "statistical" resemblance with the simulations 
more apparent; simulations are presented with the ver- 
tical and horizontal directions parallel to the computa- 
tional grid, although lattice anisotropy is not appreciable 
here. 

Patterns are shown at two different times: t = 1.32 
[Figs. Eta,c,e)] and t = 1.74 [Figs. 0>,d,f)], for all vis- 
cosity contrasts. Rather than comparing these simula- 
tions with experiments at equal dimensionless times, we 




FIG. 8: Pattern evolution for a random initial condition (see 
text) and B = 1.03X 10~ 3 , e = 0.005, dx = e, e = 1 and A = 1 
(a,b), A = 0.5 (c,d) or A = (e,f). Snapshots (3Ro X 3Ro) 
shown at t = 1.32 (left) and t = 1.74 (right). 



present all the patterns when their envelope has roughly 
attained twice [Figs. IHa,c) and[8Ja,c,e)] and three times 
[Figs. ^b,d) and|Sfb,d,f)] the initial circle radius. Direct 
comparison of the time scales is left for Sec. [V] 

Let us first compare patterns at the earlier stage [Fig. 
Ufa) with Fig. Ela), and Fig. □» with Fig. Q^b)]. 
The similarity between experiments and simulations is 
remarkable, especially if one takes into account that the 
initial conditions were only "statistically" similar (see 
above). Furthermore, the common features (overall mor- 
phology, filament accomplished or frustrated branching, 
droplet size, the already mentioned incoming fingers) are 
well reproducible in other simulation runs and experi- 
ments. This all indicates that the excellent agreement in 
the typical morphologies between experiments and sim- 
ulations is not fortuitous. This is noteworthy, since a 
change in the wetting conditions for the high viscosity 
contrast experiments, for instance, completely changes 
the observed morphology pcj . 

^From this earlier stage it would seem that pinch-off 
is not significantly more present for the lower (A = 0.5) 
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than higher (^4 = 1) viscosity contrasts. However, sim- 
ulation offers us the possibility to go to the ideal limit 
A = 0, in which both fluids have strictly equal viscosi- 
ties. Fig. [SJe) shows that most droplets have pinched-off 
by the same time for A = 0. Comparing Figs. |HJa,c,e) 
between them, it is clear that pinching arises as the vis- 
cosity contrast is decreased. 

Another possibility is to go to later times. The ex- 
periments still show no pinch-off for A = 1 [Fig. QJb)], 
but most filaments have emitted at least one droplet for 
A ~ 0.5 [Fig. IHd)]. The simulations would seem less 
conclusive, since they display some pinch-off events for 
A = 1 [Fig. Hfb)]. However, these are most likely to be 
spurious: They occur for very narrow filaments, whose 
width is comparable to the interface thickness e, and we 
know from earlier work psl l29l | that overlapping diffuse 
interfaces attract each other. Most importantly, pinching 
is inhibited as e is decreased (the same run with a larger 
value of e displays more pinching). 

Simulations for low viscosity contrasts [A = 0.5, Fig. 
|Sfd) and A — 0, Fig. Eff)], show much more pinching 
(the lower A, the more). The first pinch-off at the tip 
of a filament does not significantly change with e, very 
much like we observed for the dumbbell-shaped pattern 
in Sec. IIV Bl Some pinch-off events disappear, while 
other appear; some are reentrant in e (they disappear 
for a middle value of e but reappear when e is further 
decreased), . We hence conclude that these first pinch-off 
events are not spurious. In any case, longer and thinner 
filaments are achieved as e is decreased. 

Pearling (successive pinch-off events shortening a fil- 
ament from its outer end) is increased when decreasing 
e, even though it is less present in the experiments [Fig. 
H^d)]. It is interesting to note the absence of pinching 
at the base of the filaments, near the rotation axis, for 
A = [FiglBIf)], as opposite to the A = 0.5 and A = 1 
runs [Figs |Hl[b,d), respectively]. Since we saw pinching 
to be spurious for the A = 1 run, this suggests that such 
events at the base of a filament are also spurious. Fur- 
thermore, they are not observed in the experiments. 

We conclude that simulations reproduce very accu- 
rately the experimental morphologies before pinching, 
which strongly supports the idea that the dynamics can 
be considered to be the standard Hele-Shaw for that pur- 
pose. They also indicate that these bidimensional dy- 
namics do lead to finite-time pinching for low, but not 
high viscosity contrasts. The crossover occurs somewhere 
between A = 0.5 and A = 1. 



E. Filament stretching and droplet dynamics 

We have postponed the study of the filament stretch- 
ing and droplet scalings to present it in a unified way, 
for all the cases (n = 2, 4 and 18 filaments) considered 
above. Let us first summarize what we can expect from 
the theory in Sec. IIIII Droplets at the tip of a filament 
should run away roughly exponentially in time, with a 



rate comprised between K ' = 1 (for an isolated droplet) 
and K = 2/(1 + A) (the filament thinning rate). The 
lower bound K' corresponds to the idealized case where 
the filament does not perturb the droplet shape nor the 
flow created by it, so that the droplet still obeys the exact 
solution for an isolated and circular one. This assump- 
tion is more realistic for high viscosity contrasts, since we 
have seen that in the A — > 1 limit the homogeneous flow 
inside an isolated circular droplet holds for any closed in- 
terface, regardless of its shape and whether it is isolated 
or not. Note, however, that a droplet at the tip of a fil- 
ament is not completely closed, so that a deviation from 
K 1 is still possible. 

The upper bound K is based on another idealization, 
where no fluid is fed into the droplet, so that all the 
fluid expelled by the filament as it narrows is used to 
extend it. Therefore, the radial coordinate of the neck of 
the filament just before the droplet scales with the same 
exponent K as the filament width. If the shape of the 
droplet is stationary, the radial coordinate of the droplet 
center of mass is then just that of the neck of the filament 
plus a constant. 

The purpose of this section is to determine which is the 
actual evolution of the radial coordinate of the droplets 
center of mass. Our experiments have not clearly an- 
swered this question: Droplets are seen to (i) run away 
with roughly the same rate whether connected to a fila- 
ment or isolated, but (ii) this observed rate seems to be 
K (filament thinning) rather than K' (isolated droplet). 
The problem is that (ii) is inconsistent with the exact 
isolated droplet solution. 

Fig. [5] shows the evolution of the radial coordinate of 
a droplet center of mass for the dumbbell-shaped pattern 
in Fig. El(n = 2, dashed lines) and some of the furthest 
droplets in Fig. |H1 (n = 18, solid lines), compared to 
that of an isolated circular droplet (n = 1, dotted lines). 
Bold (thinner) lines correspond to A = 1 (0). Curves 
are shown till the end of the respective runs, except for 
the isolated droplet, which is shown only until it begins 
to "feel" the finite simulation box. Note that for n > 
2 curves make little sense for R < Rq, and this part 
has been removed for the n = 18 runs, since the code 
actually detected a different droplet at those early stages. 
More details on how the curves were obtained can be 
found below. Although the range of R/Rq monitored 
decreases with increasing n, the final aspect ratio of the 
pattern (droplet radius to length of filament for n > 2, 
or droplet radius to radial coordinate of center of mass 
of the isolated droplet) attained is similar. 

Clearly, all droplets scale with roughly the same rate 
to, although runs with A — 1 show a good linearity, 
whereas those for A — do not. More precisely, for 
A = 1 (0), one measures (in the straightest segment of 
the curves) to = 0.8 (0.7) for n = 2, to = 0.82 (0.75) 
for n = 18, and to = 0.88 (0.83) for n = 1. Runs for 
n = 4, 6, 8, 11 and 12 fingers, give similar results, also 
for intermediate values of the viscosity contrast (A = 0.5, 
0.8). 
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FIG. 9: Radial position (in log scale) of the center of mass 
of various droplets vs. time, for A = 1 (bold) and A — 
(thinner curves). Solid, dashed and dotted lines correspond 
to the droplets indicated with an "x" in Figs. |8{a,b,e,f), those 
in Fig. and the isolated droplet in Fig. 1101 bottom-right, 
respectively. Curves are shown till the end of the respective 
runs, except for the isolated droplet, which is shown only until 
it begins to "feel" the finite simulation box. Inset: Linear 
regime of the isolated droplet, continued in the main plot after 
translating the (dotted) curves to earlier times for comparison 
with the others. The inset preserves slopes with respect to the 
main plot. 



Since all droplets scale with a same rate, roughly in- 
dependent of A and close to K' = 1, it is clear that 
the droplet scaling is dominated by that of an isolated 
droplet. This rate is found to be insensitive to pinch-off, 
and it also holds for the experimental isolated droplet in 
Fig. [5] (dotted lines), which further supports this idea, 
and confirms this part of our experimental conclusions 
(i). However, the fact that droplets scale roughly with 
K' = 1 for any viscosity contrast A in our simulations 
(to ~ 0.8) and certainly not with K = 2/(1 + A) confirms 
that there is some problem with the statement (ii) that 
the experimental growth rate be K = 2/(1 + A), as al- 
ready expected from the exact isolated droplet solution. 
We address this issue in the following section. 

The mismatch between the droplet scaling (m ~ 0.8) 
and the filament thinning K = 2/(1 + A) accounts for 
the fluid fed into the droplet. This is relatively small 
for A = 1 (K = 1), but not for A = (K = 2). This 
explains the experimental and numerical observation that 
droplets grow faster and larger as A is decreased. The 
bad linearity of the ^4 = runs is also probably due 
this large mismatch, although we recall that the filament 
thinning itself displays a less constant exponent for A = 
than A = l. 

The center of mass of the droplets for n = 2, 18 in Fig. 
ED was computed subtracting the final droplet radius to 
the droplet tip position, which is not rigorous, because 
the droplet radius changes during the evolution. How- 
ever, if one subtracts, say, the final distance from the 
droplet tip to the neck of the filament before the droplet 
for comparison, rates increase by 0.1 only (not shown). 
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FIG. 10: Stream lines induced for A — by a dumbbell- 
shaped pattern as that in Fig. [S] (top; a quarter of the in- 
terface shown, B = 0.09, e = 0.0067, dx = e/2, e = 0.5) or 
by an off-center circular droplet (bottom right; a half shown, 
B = 0.32, e = 0.02, dx = e, e = 0.1, dt = 0.25ed:r 2 ). Bottom 
left: Blow-up (y enlarged 30 x more than x coordinate) of the 
forming neck in the top pattern. Interface (thicker lines) and 
contour plot of the stream function (thinner lines) at equally 
spaced levels (different in each graph); x = y — is the rota- 
tion axis. 



We have also extracted the neck position (not shown) of 
the same droplets (marked with an "x") in Fig. [3J and 
measured an exponent up to to = 0.9 (for A — 0.5). We 
also confirmed that a droplet which pinches only at the 
end of the run (marked with "+") scales with the same 
rate than one of the first to pinch (marked with "x" ) . 

In the end, the main effect of the filament on the scal- 
ing of the center of mass of the droplet seems to come 
from the injection of inner fluid into it, rather than per- 
turbing the droplet environment. This should mean that 
the flow of the inner fluid is much more affected than 
that of the outer by the presence of the filament. The 
latter is illustrated for the most sensitive case (A = 0) in 
Fig. 1101 where stream lines are superimposed to the in- 
terface for a quarter of a dumbbell-shaped pattern (top), 
as that in Fig. and for half an isolated droplet (bot- 
tom right): The most remarkable difference is the higher 
density of stream lines (higher velocity) at the entrance 
of the droplet and inside it, which breaks the flow uni- 
formity within the droplet; the outer flow is qualitatively 
similar near the droplet (note that stream lines are not 
shown at the same levels for a maximum resolution of 
the flow details in both cases, so that the modulus of the 
velocity cannot be compared). 

When a neck is formed (blowup, bottom left, in Fig. 
llOp . the flow is affected also outside the droplet, but 
mostly near the neck. Because pinching then also occurs 
in a time scale sensibly shorter than that of the filament 
stretching, the effect is not perceived in the latter. The 
A = droplets tracked in Fig. for instance, pinch dur- 
ing t € (0.9, 0.96) for n = 18 and t € (2.3, 2.4) for n = 2, 
but they curve down well after or before it, respectively. 

One could also wonder why the exponents measured 
for the isolated droplet run are not strictly to = K' = 1, 
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but this is certainly a finite- interface-thickness correction 
due to the less refined value used (e = 0.02). Our purpose 
here was merely to check that the exact solution for the 
isolated droplet has some finite basin of attraction, as 
we find it to be the case: The droplet stays circular and 
runs away exponentially with a very steady rate from 
the linear (q = 0.005, inset of Fig. |5J) to deep into the 
nonlinear (main plot) regimes. The droplet depicted in 
Fig. bottom right corresponds to the end of the linear 
(and beginning of the nonlinear) regime. The type of flow 
predicted by the exact solution [Eqs. {HJ], is apparent 
here. 



TIME SCALE 



In Sec. IIV El we saw that the statement that the exper- 
imental growth rate be if = 2/(1 + A) and not K ' = 1 is 
inconsistent not only with the exact isolated droplet so- 
lution, but also with numerical simulations of the Hele- 
Shaw equations. This statement relies only on the ob- 
served A = 0.45 rates, since for A = 1K = K' = 1 (in 
good agreement with the measured values). This indi- 
cates that our experiments for low A do not follow the 
standard Hele-Shaw dynamics. 

A three-dimensional study of the proper effective 
boundary conditions to apply on an idealized two- 
dimensional interface, as performed for air displacing a 
liquid [U, is lacking for the case of two liquids. This 
low-A case might be different, since a second viscosity 
is involved when a layer of vaseline oil (outer fluid) wets 
the plates also inside the silicone oil (inner fluid) domain. 
Note that, for A = 1, such a wetting layer (in this case, 
of silicone oil) is to be found in the air. However, differ- 
ent boundary conditions should in general not produce 
experimental patterns so similar to the outcome of sim- 
ulations using standard boundary conditions as those in 
Sec. IIVDI Especially, they should not necessarily pre- 
serve the pinch-off for low A. 

The inferred dimensionless rates for the experimental 
droplet runaway rely on the measured rates and on the 
time scale t*. A modification of the latter can hence al- 
ready account for the observed discrepancy between the- 
ory and simulations. To further realize this, it is instruc- 
tive to compare again experimental (Fig. ^) and simula- 
tion (Fig. |SJ) patterns: The similarity in the morpholo- 
gies is obvious, but time scales are not straightforward 
to compare, since the degree of deviation from a perfect 
circle in the experimental initial condition is unknown. 
This results in an effective "latency time" , during which 
no significant departure from a circle is observed. 

However, we can compare the time elapsed for the pat- 
tern envelope to grow from a radius of 2i?o [Figs. QIa,c) 
and Efocje)] to 3i? [Figs. QJb.d) andEtb,d,f)]. This is 
At =6s(36s) for A=l(0.45) in the experiments. With the 
measured values of the physical parameters, t* in Eq. ^ 
is t* = 15.5s(134s), which results in a dimensionless time 
interval At/t* = 0.39(0.27) for A=l(0.45). This is to 



be compared with At/t* = 0.42 VA in the simulations. 
The mismatch between experiments and simulations is 
not significant (taking into account the viscosity uncer- 
tainty and the subjectivity in the snapshots to compare) 
for A = 1, but it is for A — 0.45. This is just a (sta- 
tistical) confirmation of the droplet scaling, since it is 
this scaling what determines the growth of the pattern 
envelope, but makes it clear that the main difference be- 
tween experiments and simulations is the time scale for 
A = 0.45. 

A simple explanation is that the vaseline wetting layer 
has a more obvious effect than modifying the effective 
two-dimensional boundary conditions on the interface: it 
changes the expression for Darcy's law [Eq. (fSJ] . The 
two-dimensional velocity in Darcy's law results from an 
average over the cell gap of the three-dimensional veloc- 
ity field. We simply point out that one should average 
over the fraction of the gap actually filled with silicone 
oil when computing the averaged silicone oil velocity. As- 
suming stick boundary conditions for the vaseline on the 
glass plates, and no-slip between vaseline and silicone, 
we average the silicone velocity in a region closer to its 
maximum (at the midplane between the two glass plates) 
than in the absence of the vaseline wetting layer. This 
results in a higher velocity, further enhanced by the fact 
that vaseline is less viscous. More precisely, we find 



b 2 [l + Ab 2 ] 

12/Xi n 



where 



Ab 2 = ^ + I 8 + 6 



A*out 



rw\ 2 



(16) 



(17) 



with w the thickness of the vaseline wetting layer at 
each glass plate. The effect of this on the time scale 
of the dynamics t* in Eq. is to replace /j,- ln there by 
/ii n /[l + Ab 2 } (wout stays unchanged). For a wetting layer 
of thickness w 10% of the total cell gap b, our A = 0.45 
time scale should decrease by a factor 0.78. This re- 
sets the rate at which droplets run away to m = 0.97- 
1.17, in reasonable agreement with theory and simula- 
tions. We have not measured the thickness of the wetting 
layer w, and possibly other effects like the effective two- 
dimensional boundary conditions on the interface might 
play a role, but this gives a plausible explanation of our 
observations. 



VI. LUBRICATION THEORY 

The purpose of this section is to derive the time evolu- 
tion of a gently curved filament in a more systematic way. 
In particular, we would like to account for the striking 
difference between the filament thinning of high and low 
viscosity contrast dynamics. 

To do so, we perform a lubrication approximation for 
the interface. This assumes the interface height h(x,t) 
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in Fig. |3 to vary in a scale Ah much smaller than the 
scale of horizontal variations £, and expands formally all 
quantities in powers of e = Ah/£. We then find an evo- 
lution equation for h up to first order in e. The idea is 
to see whether the viscosity contrast can affect the rela- 
tive importance of the various terms and to what order 
it enters. 

We begin by rewriting the free boundary problem of 
Eqs. <|14(l into an exact evolution equation for h: Using 
the vortex-sheet formalism of Refs. |26|,|32], we have 



d t h = Uy - U x h x , 



(18) 



where U(x, t) = U x x + Uyi) is the mean fluid velocity 
at the interface, which we assume to vary only in the x 
direction: 



an arbitrary way. We pursue the calculation and check 
whether the result is independent of the shape of rj. This 
cut-off is necessary, as opposed to the channel case, be- 
cause the zero order of the vorticity in a rotating cell 
is neither zero nor a constant but, Y^ ) oc x. The other 
major difference between the weakly nonlinear expansion 
performed in and this lubrication approximation is 
the appearance of the second integral in (|19|l . Its ex- 
pansion is not trivial as pointed out in [T^j : there it was 
solved using limiting procedures, the Plemelj formula, 
and delta function representations. 

Substituting these results into Ea. (|21f) . undoing only 
the scaling of the velocity, and dropping the cut-off func- 
tion rjlx), which is not necessary for finite or asymptoti- 
cally straight infinite filaments, we finally obtain 



U(t — 1 P f + °° (h(x')-h(x),x~x') ~V n , / 

1 p r+°° {-h(x')-h(x),x-x') 
2tt J-oo (x-x') 2 + (h(x')+h(x)y 2 1 V x > 



(19) 



where P is the Cauchy principal part and T = T(l + h 2 ,) 
(7defined in Sec. IIV A (I becomes 

T/2 = Bd x K -AU-(x + yd x h) +x + hd x h. (20) 

We now scale h with Ah, x with £, and the interface 
velocity U with Vo = (B/i 2 )+£. We expand any quantity 
Q as Q = Q (0) + eQ (1) + e 2 Q (2) , so that the evolution 
equation for h up to O(e) becomes 



dth = U 



(1) -B x wf +e(llf ) -d x hVf\ (21) 



where we have anticipated that = 0. U and T are 



indeed expanded along the same lines as in Refs. |14L I3S 
to further find 



U 
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UP = + H[d x (hT^)] (22) 



f ] = -lldx(hT^) + hd x r^} 



(23) 



U 



f ] = -l[d x (hT^) + hd x r^] _ h H[d 2 x (hT^)}, (24) 



T (o) 
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r/(x)x 



(25) 



T (i) = 



- B M-J^r H[d x (h?W)}, (26) 



(1 + A)V a L 2 



1 + A 



where H(x) = n~ l P J+~ f(x')dx'/(x-x') is the Hilbert 
transform of f(x), and r](x) is an arbitrary cut-off func- 
tion in the centrifugal force to render H finite, which is 1 
up to a certain distance xf and then decreases to zero in 



^d t h 



-d x {xh) - eBd x {hdlh) 
-e^§d x (hH[d x (xh)]), 



(27) 



We see that this is a simple continuity equation for h, 
namely that already used in Sec. IIII Al d t h = —d x (hv x ) 
[Eq. (0)], but now with an inner fluid velocity along the 
x direction 



1 + A 



sBdlh 



e l —^H[d x {xh)\ 



(28) 



which is no longer just v x = Kx as used in Sec. IIII Bl but 
takes into account all the terms entering the tangential 
velocity jump there [Eq. |J3J]. 

We now present a general solution when the order e in 
Ea. (|77|) is neglected. Any given initial condition ho(x) 
evolves following a simple scaling h(x, t) = Lho(Lx) [Eq. 
©] with L = e~ Kt and K = 2 /(1 + A). This agrees with 
the scaling found in Sec. IIII Al for an initially straight in- 
terface. Here, we find that the result is more general and, 
as a consequence, any interface approaches zero at least 
at infinite time. In this way, rotation produces pinch- 
off at infinite time independently of the initial interface. 
This is to be compared with the influence of a gravity jet 
which only produces a shift in the moving frame. 

This leading behavior smooths out the gradients both 
by lowering height differences [Eq. JSJ] and by stretching 
neighboring points apart [Eq. ©]. One would expect 
straight filaments to be stable, and the solution to be- 
come increasingly accurate in time. This smoothing ef- 
fect of the centrifugal force through the zeroth order term 
v x = Kx presumably competes with the higher order 
terms: The term in dimensionless surface tension B, for 
instance, is known to lead to finite-time pinch-off in the 
absence of any other force and for an inviscid outer fluid 
(A = 1), for certain initial conditions [lCl llll Il4|. The 
analysis of the interplay between rotation and surface 
tension and their competing effects deserves an in-depth 
separate analysis that will be carried out elsewhere. 

However, we can address in this framework the system- 
atic observation of pinch-off events in experiments and 
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simulations for low viscosity contrasts A = 0-0.5 and the 
lack of them for high viscosity contrast. This fact cannot 
be explained by the two first terms in the r.h.s of (|27|1 . 
If such singularities and differences are contained in the 
Hele-Shaw dynamics as the simulations seem to indicate, 
they should be linked to the remaining non local term 
of order e. Although this term does indeed play a role 
for A < 1, we have seen that it enters at the same or- 
der than the surface tension term. The question is then 
how this, in principle, higher order term can affect the 
general asymptotic scaling of filaments which typically 
become straight for any viscosity contrast. 

The answer is in the non-local nature of the new term. 
Hence, a locally straight filament is still influenced by a 
curvature elsewhere. Real filaments are finite, and this 
means that the upper and lower interfaces necessarily 
meet in a highly curved region (in practice, a droplet at 
the tip), where the e expansion breaks down. Because 
this will be felt by an elsewhere straight filament, the e 
expansion is much more delicate for low viscosity con- 
trast, and significant deviations from the scaling of Eqs. 
© or (J5J can be expected. The influence of a curved re- 
gion decays with distance, so that these deviations should 
decay for the central part of a filament as it grows long 
enough for its ends to have little effect on it. 

Although it is difficult to address the precise effect of 
the non-local term analytically, we can gain some in- 
sight by considering the simplest possible situation: a 
perfectly straight filament of time-dependent arbitrary 
width. More precisely, we consider only an upper straight 
interface given by y — h(t) in Fig. [3] and the region 
y > 0, x > 0, with no-flux boundary conditions on the 
axis y = and x — 0, exactly as in our simulations in 
Sec. IIVBI Proposing a stream function ip (bi)linear in x 
and y to solve the full Hele-Shaw equations [Eqs. I|14|) ]. 
we find 

ipin = K- m xy, (29a) 
i/w = h(K in - K out )x + K QUt xy, (29b) 
h = h(t = 0)e- K - 4 , (29c) 

K in = K+ \ / A out , (29d) 
1 + A 

where we recall that A = 2/(1 + A). 

For A = 1, we recover K ln = A because the inner 
fluid and interface dynamics are completely specified and 
decoupled from the outer flow. For A < 1, however, the 
whole solution depends on the outer flow, and hence on 
actual external boundary conditions. For a semiinfinite 
filament and no external breaking of the translational 
symmetry, A out = 0, and the scaling solution is again 
valid and exact. 

For a finite filament, though, there will necessarily be 
a recirculation (not to speak of the flow created by a 
droplet): Stream lines going out of the filament at its tip 
will come back to different points of the interface, in a 
loop of the size of the filament length R itself. The outer 
flow solution proposed will no longer be valid, but, for 



a long enough filament, we can assume that the solution 
still holds for x -C 1. Given the fact that recirculation 
requires a positive hyperbolic profile close to the origin 
(see Fig. llUfl . A' out must be positive and therefore K m > 
A, slowly approaching A as the filament grows longer. 
This is consistent with the filament thinning exponent 
for A = observed in Fig. 

One could speculate that this effect becomes stronger 
as one approaches the filament tip (although the flow 
also looks different), making the local, effective thinning 
(and its rate) increase. This would lead to the forma- 
tion of a more shallow region (the neck) and thus initiate 
the pinch-off. The forming neck then happens to receive 
a much higher density of incoming stream lines that its 
vicinity (see Fig. ^3 blowup at the bottom left), pushing 
it further down. Surface tension enters the problem as 
the filament is locally curved, so that it is not possible 
to tell whether the pinch-off itself is due to the non-local 
term or to surface tension, or yet to a combination of the 
two unless a careful analysis is performed. It is neverthe- 
less clear that the non-local term does provide the initial 
departure from the scaling solution needed to overcome 
its smoothing effect. 



VII. DISCUSSION AND CONCLUSIONS 

We have studied viscous fingering in a rotating Hele- 
Shaw cell from the point of view of the dynamical ap- 
proach to pinch-off singularities. We have found that this 
is linked to a global scaling of the pattern: Radial fila- 
ments of denser inner fluid narrow and stretch; they push 
away the droplet at their tip, and incoming fingers of less 
dense, outer fluid approach the rotation axis to compen- 
sate for the flux expelled through the filaments. For n- 
fold symmetric patterns, the system can approach cither 
a steady shape in infinite time or a finite-time pinch-off 
singularity. The former consists of n convex arches end- 
ing in straight, formally semiinfinite radial filaments of 
zero width with a droplet at their tip (at infinity) con- 
necting each pair of arches and running away with infinite 
velocity. 

Thus, we see that rotation enforces pinch-off as ex- 
pected, at least in an infinite time. More precisely, at 
the lowest order in thickness variations of a lubrication 
approximation for the filaments, they narrow exponen- 
tially in time. Incoming fingers slow down in the same 
way, whereas the filaments stretch and the droplets at 
their tips are centrifuged away also exponentially. 

In the high viscosity contrast limit, A — > 1, all these 
exponential behaviors have the same time constant, l/t*, 
where t* is the characteristic time scale of the problem, 
and we have been able to observe all of them and confirm 
this value by phase- field simulations and experiments. 
For low viscosity contrasts, the filament thinning has a 
time constant of A = 2/[(l + A)t*] in absolute value, 
whereas the droplet dynamics remain close to A' = l/t* 
as for high viscosity contrasts. The droplet scaling is 
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thus universal in a first approximation, i.e., only weakly 
dependent on any dimensionless control parameter. Fur- 
thermore, we have found it not to be much affected by 
whether the droplet is isolated or connected to a filament. 
We have derived an exact time-dependent solution (pre- 
viously found in the A — ► 1 limit |2J|) for the case of an 
off-center isolated circular droplet which shows that this 
universality is rigorous for isolated circular droplets. 

In a second approximation, droplets at the tip of the 
filaments are actually found to scale with a time constant 
slightly below the theoretical value 1/t*, so that this is 
always below that of the filament thinning 2/[(l + A)t*]. 
Some fraction of the fluid expelled by the filament as it 
narrows is hence fed into the droplet. Because the time 
constant for the thinning increases above as A is 
decreased, the lower A the more flux the droplets receive, 
as observed in both experiments and simulations. 

For experiments with two liquids (low A), comparison 
with the exact solution for an isolated circular droplet 
and extensive simulations have shown the time scale of 
the dynamics t* to be lower than its characteristic value 
t* = 12(/ij n + /i out )/(& 2 £l 2 Ap). We explain this by an ef- 
fective lower viscosity of the inner fluid /i; n , which comes 
from Darcy's law corrected to take into account the pres- 
ence of a wetting layer of the outer liquid. This results 
in an increased inner liquid mobility. 

Although rotation causes infinite-time pinch-off at the 
lowest order in the variations of the filament thickness, 
both the thickness decay and the stretching apart of 
neighboring points smooth out thickness variations. Ro- 
tation also enters the other next order term: a non-local 
term proportional to 1 — A. Because this term is non- 
local and the interface is necessarily highly curved in the 
droplets at the end of the filaments, it can overcome the 
stabilizing effect of the lowest-order contribution of rota- 



tion, and thus initiate finite-time pinch-off. The role of 
surface tension and the interplay with the local and non- 
local terms associated with rotation remain interesting 
open questions. 

In any case, this non-local term can account for the dif- 
ferences observed in both experiments and simulations, 
where finite-time pinch-off singularities appear system- 
atically for low values of A but not for high values of 
A. Three-dimensional effects like wetting do not seem to 
account for the initiation of pinch-off, since simulations 
are free of them and nevertheless reproduce well the ex- 
perimental morphologies and the pinch-off phenomenon 
itself. Corrections due to the non-local term are also the 
main reason for obtaining a less clear exponential behav- 
ior for low viscosity contrasts. 

Note that, for the case of a channel under gravity or 
pressure-driven injection, this non-local term still arises, 
although it has another form. Therefore, it may also 
explain the enhancement of finite-time pinch-off observed 
for low A in these well known cases. 
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